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Abstract 

To  provide  answers  to  the  concern  as  to  how  quickly  the  temperature  of  solid  oxide  fuel  cells  (SOFCs)  for  transportation  application  will 
drop,  a  thermal  analysis  of  the  cool-down  time  of  an  SOFC  stack  during  vehicle  idle  or  stand-by  has  been  carried  out.  Because  a  large  amount 
of  thermal  energy  is  stored  in  high-temperature  SOFC  stacks,  it  is  important  to  select  suitable  thermal  insulations  to  reduce  heat  loss.  Three 
typical  kinds  of  thermal  insulating  materials  have  been  selected  in  the  present  calculations.  The  results  indicate  that  a  high-performance, 
vacuum-multifoil  thermal  insulation  can  be  applied  to  significantly  reduce  heat  loss  and  to  maintain  temperature  uniformity  across  a  cell  stack. 
Consequently,  the  cool-down  time  from  1000  to  800  °C  is  extended  from  2  h  (with  a  5  cm  thick  conventional  material)  to  about  31  h  (with 
a  1  cm  thick  high-performance  material). 
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1.  Introduction 

Solid  oxide  fuel  cells  (SOFCs)  have  many  attractive  fea¬ 
tures  [1,2]  (e.g.  high  energy  efficiency,  low  emission,  rapid 
refueling,  etc.),  and  are  thus  under  consideration  as  an  alter¬ 
native  power  source  for  transportation  application  [3],  Cur¬ 
rent  SOFCs  must  be  operated  at  temperatures  around  1000 
°C,  in  order  to  reduce  the  resistance  of  yttria  stabilized  zir- 
conia  (YSZ)  electrolytes.  Many  efforts  [4]  have  been 
focused  on  reducing  SOFC  operation  temperature  from  1000 
°C  to  as  low  as  500  °C,  by  developing  thinner  YSZ  layer 
fabrication  technology  [5]  or  by  developing  electrolytes  with 
lower  resistance  [  6] .  However,  even  for  the  lower  operating 
temperatures  of  about  500  °C,  thermal  management  plays  an 
important  rote  in  the  successful  operation  of  SOFCs. 

In  general,  a  thermal  management  system  is  responsible 
for  raising  the  temperature  of  an  SOFC  stack  to  its  operating 
temperature  before  system  start-up,  and  for  maintaining  the 
operational  temperature  and  temperature  uniformity  across  a 
stack  during  cell  operation  (e.g.  if  excess  heat  is  generated 
during  ceil  operation,  cooling  is  required,  also,  thermal  insu¬ 
lation  is  needed  during  vehicle  idling  or  stand-by).  In  typical 
SOFCs,  the  thicknesses  of  the  cathode  and  the  anode  are  of 
the  magnitude  of  about  100  to  200  (un,  the  electrolyte  thick- 
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ness  is  of  a  few  tens  of  micrometers,  and  fuel  and  air  gases 
flow  through  the  fuel  and  air  channels  in  the  interconnect 
layers.  Therefore,  the  warm-up  time  may  be  kept  short 
(because  thermal  energy  can  be  brought  into  the  cell  by  hot 
gases,  and  the  heat-transferpath  is  short  within  each  individ¬ 
ual  cell ) ,  if  the  all-ceramic  structure  of  SOFCs  can  withstand 
large  thermal  stresses.  The  maintenance  of  operating  temper¬ 
ature  and  temperature  uniformity  during  cell  operation  will 
need  judicious  design  of  the  cell  stack  including  cell  config¬ 
uration  ( cross-flow,  co-flow  and  counter-flow) .  For  the  trans¬ 
portation  application  of  SOFCs,  one  important  concern  is 
how  quickly  the  temperatures  in  a  stack  will  drop  below  the 
operating  temperature  at  open-circuit  conditions,  which  is 
related  to  the  requirements  for  thermal  insulation.  In  this 
report,  a  simple  thermal  analysis  (for  a  fuel  cell  stack)  of 
cool-down  time  and  suitable  thermal  insulation  is  presented. 
In  addition  to  the  fuel  cell  stack,  a  SOFC  power  source  would 
also  contain  a  fuel  processor,  an  air  preheater,  and  some  other 
components.  Therefore,  more  rigorous  analyses  will  have  to 
be  done  after  a  whole  system  has  been  engineered. 

2.  Planar  SOFC  configuration  and  the  mathematical 
model 

A  tubular  design  of  SOFCs  is  now  at  the  most  advanced 
stage  of  development  among  different  designs.  However,  pla- 
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Fig.  1 .  Schematic  setup  of  a  planar  SOFC  stack  and  unit  cell  dimensions. 


nar  SOFCs  are  considered  to  be  more  compact  and  of  higher 
power  density  in  comparison  with  the  tubular  SOFCs.  There¬ 
fore,  a  planar  SOFC  configuration  is  adopted  in  the  present 
modelling.  A  planar  SOFC  stack  [7]  is  shown  in  Fig.  1.  A 
stack  includes  many  cell  layers  and  each  layer  consists  of 
many  unit  cells.  A  cell  is  composed  of  three  layers:  Ni0  35 
(Zr02)o.65  cermet  anode  ( fuel  side  electrode,  porosity  0.4  to 
0.6),  (Zr02)o.9o  (Y2Oj)01o  electrolyte  (porosity  <0.06), 
and  Lao.^Sro.^MnOj  cathode  (air  side  electrode,  porosity 
0.3  to  0.5).  The  ceramic  interconnects,  Lao.K4Sr0|6Cr03 
( porosity  <  0.05 )  act  as  bipolar  current  collectors  ( if  a  high- 
temperature  alloy  or  a  conducting  oxide-coated  metal  is  used 
as  bipolar  plate,  the  present  thermal  calculations  should  be 
conducted  with  a  suitable  adjustment  of  thermophysical  par¬ 
ameters).  Fig.  1  also  shows  the  unit  cell  dimensions  [8].  To 
obtain  high  power  and  to  simplify  the  design  for  bringing  air 
and  fuel  into  SOFCs,  stacks  with  large  cross-sectional  areas 
are  preferable  (provided  that  additional  problems  will  not  be 
caused  by  temperature  and  current  distribution).  However, 
only  cell  sizes  of  up  to  400  cm2  have  been  fabricated  to  date 
[9] .  Therefore,  in  this  calculation,  the  cross-sectional  area  of 
the  stack  is  assumed  to  be  20  cmX20  cm  (i.e.  Ly=Lt- 20 
cm,  see  Fig.  1).  With  current  SOFC  technology,  electrode 
power  densities  below  0.25  W  / cm2  have  been  achieved  1 10] . 
If  a  power  output  of  20  kW  is  designed  for  an  SOFC  system, 
a  minimum  total  electrode  area  of  8  m2  is  required,  that  is, 
200  cells  with  electrode  areaof  20  cm  X  20  cm.  Consequently, 
the  thickness  Lx  of  the  stack  is  86  cm  ( thermal  insulation  not 
included). 

An  SOFC  stack  has  heterogeneous  thermal  properties.  It 
is  reported  [11]  that  if  this  characteristic  is  taken  into  con¬ 


sideration,  the  time  to  compute  temperature  distributions  will 
be  too  long  for  stack  s  of  practical  size,  which  is  similar  to  our 
observation  in  thermal  modelling  of  lithium/polymer  batter¬ 
ies  [12].  A  way  of  circumventing  this  difficulty  is  to  treat 
the  stack  as  having  effective,  anisotropic  thermal  properties. 
This  treatment  is  a  good  approximation  for  a  stack  consisting 
of  many  cell  units  and  many  cell  layers.  In  this  case,  the 
transient  heat-transfer  equation  for  an  SOFC  stack  at  open- 
circuit  conditions  can  be  represented  as: 


a  t  a2r  a2r  a 2r 


where  kx,  ky,  k.  are  the  effective  thermal  conductivities  in 
different  directions,  p  and  Cp  are  the  average  stack  density 
and  specific  heat,  respectively. 

Heat  loss  is  caused  by  natural  air  flow  and  thermal  radiation 
on  the  external  surface  of  a  stack.  To  significantly  reduce  heat 
loss  to  the  surroundings,  high-performance  thermal  insulation 
materials  may  be  needed  so  that  the  temperature  at  the  exter¬ 
nal  surface  of  the  insulation  case  is  close  to  that  of  the  ambient 
air.  It  is  assumed  that  an  SOFC  stack  is  enclosed  in  an  insu¬ 
lation  case  and  that  the  temperature  gradients  within  the  insu¬ 
lation  case  will  be  at  a  quasi-steady  state.  Therefore,  the 
boundary  condition  for  thermal  insulation  established  for  lith¬ 
ium/polymer  batteries  [  13]  can  be  used,  that  is: 


-k^\n,0M  =  (hn+hM\+'I^)-'(Tn-T*)  (2) 


where  fir=4  t<xT£  [14]  is  a  radiation  heat-transfer  coeffi¬ 
cient.  hn  the  coefficient  of  heat  transfer  by  natural  air  flow  in 
different  directions,  L„  (n=x,y,z)  the  dimension  of  the  cell 
stack,  Ij„  the  thickness  of  thermal  insulation  case,  kia  the 
thermal  conductivity  of  the  insulation  case,  and  Tx  the  ambi¬ 
ent  temperature. 

The  initial  condition  is: 


T=  T0  at  /=0,  any  x,  any  y,  and  any  z  (3) 


where  7],  is  the  operational  temperature  of  an  SOFC  stack. 
The  above  equations  have  been  solved  by  the  modified 
implicit  alternating  direction  technique  [  15]  on  a  DEC  5900 
RISC  computer. 


3.  Thermo-physical  parameters 

It  appears  that  the  thermo-physical  properties  for  the  com¬ 
ponent  materials  of  SOFCs  are  not  available  in  the  open 
literature.  In  the  present  work,  the  thermo-physical  properties 
of  anode  (Ni035(ZrO2)065),  cathode  (LaaS4Sr(U6Mn03), 
electrolyte  ((Zr02)o..;o(Y203)(no),  and  interconnect 
(La,)MSro  if,Cr03)  are  estimated  in  terms  of  the  properties 
for  Ni,  ZrO:,  La203,  Sr203,  Mn203,  Mn203,  Y203,  Cr203, 
and  the  compositions  of  the  cell  components  (composition 
data  in  volume  percentage  is  used  for  p  and  k,  while  weight 
%  for  Cp ) .  The  properties  for  some  of  the  oxides  are  unknown 


Y.  Chen,  J.  W.  Evans  /Journal  of  Power  Sources  58  ( 1996)  87-91 


89 


Table  I 

Thermo-physical  properties  of  nickel,  some  oxides,  air  and  fuel 3 


Parameter 

Ni 

ZiO, 

La,0, 

Sr20, 

Mn20, 

YA 

CrA 

Air 

h2 

p(g/cm3) 

8.90 

5.60 

6.51 

4.63  hz 

4.50 

5.01 

5.21 

0.353X  I0~3 

0.0244  XlO'3 

Cp(J/(gK)) 

0.5178 

0.6149 

0.4061 

0.5366 

0.8694 

0.5621 

0.8318 

1.10 

14.71 

*(W/(cmK)) 

0.718 

0.00197 

0.03 c 

0.0215 d 

U.G34 

0.0374 

001 6'5 

0.672  X  10’J 

4.28  X IO-3 

*  Collected  from  Refs.  [  16-20) .  Those  values  are  considered  as  the  average  values  within  the  temperature  range  of  interest. 
h  Averaged  value  of  SrO  (density  4.70  g/cm3)  and  StO,  (density  4.56  g/cm3). 
c  Arbitrarily  chosen  from  the  value  range  for  ceramics. 
d  Value  for  SrO + Ti02  +  Li,0  -  ZiO,  -  Si02. 

‘  Value  for  Cr  A  +  MgO  +  Fe2Ov 

Table  2 


Thermo-physical  properties  of  anode,  cathode,  electrolyte,  interconnect  and  an  SOFC  stack 

Material 

p(g/cm3) 

Cp(J/(g  K>) 

*(W/(cmK)) 

Nlnw(Zi0>)„M  (porosity  0.50) 

3.03 

0.5951 

0.0584 

Lafl*4Sr01AMnO,  (porosity  0.40) 

3.31 

0.5726 

0.0186 

(Zr02)o.9o(YA)o.,o  (porosity  0.06) 

5.16 

0.6058 

0.0216 

Lao.a4Sr0.i„Cr03  (porosity  0.05) 

5.55 

0.5562 

0.0228 

SOFC  stack 

p(g/cm3) 

C„(J/g  K) 

UW/cmK) 

3.54 

0.5582 

*,=0.01132 

*,=0.01606 

*4=0.01444 

too,  and  are  chosen  within  a  reasonable  value  range  for  ceram¬ 
ics  as  indicated  in  Table  1 .  Table  2  lists  the  estimated  values 
for  the  thermo-physical  properties  of  cell  components  and 
the  average  density,  specific  heat,  and  the  effective  thermal 
conductivities  of  a  cell  stack.  The  latter  are  determined  in  a 
similar  way  to  calculating  the  equivalent  electric  resistances 
in  an  electric  circuit  (series/parallel  models). 


4.  Results  and  discussions 

In  terms  of  the  heat  capacity  (pCp=  1.976  J/(cm3  K)) 
and  the  size  (L,=86  cm,  Ly-L.=  20  cm)  of  the  20  kW 
SOFC  stack,  the  thermal  energy  stored  in  the  1000  °C  stack 
is  estimated  to  be  18.4  kWh  (with  respect  to  25  °C),  which 
is  of  the  same  magnitude  as  the  electric  energy  output  of  the 
20  kW  stack  during  1  h  operation.  In  other  words,  every  100 
°C  decrease  of  the  stack  temperature  will  correspond  to  1.89 
kWh  thermal  energy  loss.  Therefore,  it  is  important  to  use 
suitable  thermal  insulation  to  reduce  heat  loss  (a  significant 
amount  of  heat  is  dissipated  to  the  surroundings  in  the  form 
of  thermal  radiation,  in  addition  to  the  heat  loss  through  the 
natural  air  flow  on  the  stack  surface). 

The  operating  temperature  of  conventional  SOFCs  is  about 
1000  °C.  However,  cell  operation  may  be  started  up  at  800 
°C,  and  it  is  possible  to  raise  the  stack  temperature  from  800 
to  1000  °C  by  means  of  the  heat  generated  during  cell  oper¬ 
ation.  Therefore,  in  the  present  calculations,  the  cool-down 
time  is  defined  as  the  time  period  during  which  the  tempera¬ 
ture  at  the  center  of  a  stack  drops  from  1000  to  800  °C.  Also, 
the  cool-down  time  of  a  stack  from  1000  to  500  °C  is  calcu- 


Ipflll7  Wm-'K* 


Fig.  2.  Temperature  variation  as  a  function  of  time;  SOFC  stack:  L, = 86  n), 
L»=Z,= 20 era.  T0=  1273  K;  thermal  insulation  material:  *=0.17  W/(m 
K). 

lated,  because  many  activities  are  under  way  aimed  at  devel¬ 
oping  thin  film  SOFCs  with  lower  operating  temperatures  of 
around  500  °C. 

Cool-down  time  is  dependent  upon  the  thermal  insulation 
materials.  Conventional  thermal  insulations  are  of  thermal 
conductivities  ranging  from  0.17  to  0.03  W/(m  K).  Fig.  2 
shows  the  variation  of  the  temperatures  at  the  center  and  at 
one  of  the  eight  comers  of  a  stack  ( inside  the  insulated  case) 
as  a  function  of  time  when  a  vehicle  is  at  stand-by.  Calcula¬ 
tions  for  two  different  thicknesses  (5  and  10  cm)  of  the 
material  (fc=0. 17  W/(mK))  are  conducted.  The  cool-down 
time  for  5  cm  thick  insulation  is  about  2  h  (to  800 °C)  and  5 
h  (to  500  °C),  respectively,  while  for  10  cm  thick  insulation, 
the  cool-down  time  is  3  h  (to  800  °C)  and  8  h  (to  500  °C), 
respectively.  The  temperatures  in  the  end  cells  are  much 
lower  than  those  of  the  cells  in  the  stack  center.  Therefore, 
better  thermal  insulation  is  needed. 
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